Simulation of Many-Body Hamiltonians using Perturbation Theory with 

Bounded-Strength Interactions 



Sergey Bravyi David P. DiVincenzo t, 1 Daniel Loss * , 2 and Barbara M. Terhal ^ 

'IBM Watson Research Center, P.O. Box 218, Yorktown Heights, NY, USA 
Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4-056 Basel, Switzerland 

(Dated: March 18, 2008) 
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| Hamiltonian which relies on the Schrieffer- Wolff transformation of many-body physics. 
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INTRODUCTION 



In quantum field theory and quantum many-body physics perturbation theory is omnipresent. Perturbation theory 
is perhaps most important as a method for constructing an efficient description of the physics at low energies as 
an effective theory of the physics at higher energies. For systems which do not have a perturbative treatment, it 
can be much harder or impossible to obtain such a concise and efficient description. In an apparently independent 
development, perturbation theory has recently appeared as a tool in the area of quantum complexity theory [lj. Its 
c~| . use in this area resembles that in theoretical physics: it relates the low-energy properties of different Hamiltonian 
models. But in computer science this association is made in the reverse direction as compared with theoretical physics. 
In many-body physics, the "full" Hamiltonian (including high energy degrees of freedom) is given, and an effective 
Hamiltonian on the low-energy subspace is calculated using perturbation theory. In computational complexity theory, 
the effective Hamiltonian is given - it is a "target Hamiltonian" H talgc t chosen for some computational hardness 
property (for example, in adiabatic quantum computation Q the ground state of .fftarget encodes the execution of 
a quantum algorithm). Then a "high-energy" simulator Hamiltonian H is to be chosen such that -fftargct can be 
qq | obtained from it by perturbation theory. The main objective is to make the simulator Hamiltonian H as simple and 
\Q . realistic as possible while retaining the computational hardness of -fftargct- For example, the perturbation theory has 
' been used in Q to prove universality of quantum adiabatic computation with local two-body Hamiltonians on a 2D 
^T) . square lattice. For more recent developments see @, 0, H, 0j- 

In all applications mentioned above, one deals with a many-body Hamiltonian H — H + V where H is a simple 
operator with a known spectrum and V is a perturbation. The spectrum of H consists of two disjoint parts spanning 
the low-energy subspace P and the high-energy subspace Q. The two parts of the spectrum are separated from 
each other by an energy gap A. A perturbation expansion provides a systematic way of constructing an effective 
Hamiltonian H e g acting on the low-energy subspace P such that the ground state energy of H c g approximates the 
one of H with an arbitrarily small error. An important shortcoming of the perturbative approach is the limited range 
of parameters over which it can be rigorously justified, namely, ||V|| < A/2, see for instance [H, S|- For ||V|| A, 
we expect states from the high-energy subspace to mix strongly with states from the low-energy subspace. Thus, it 
should not be expected that a general perturbative expansion for H e g is convergent when ||V|| S> A. (One notable 
exception to this generic behavior is the case when the low-energy subspace is one-dimensional, see [9j.) 

If we require for convergence that MV"JI <C A, we force A to scale with system size n, since \ \V\ \ typically scales with 
n; this is what has been done in 0, H, [j, 0, H, ■ Hence these calculations would indicate that perturbation theory 
can only be rigorously applied to an unphysical Hamiltonian, one for which the gap A grows with the system size n. 

If the applicability of perturbation theory were thus limited, its common use in physics would be unwarranted. 
The lack of convergence of the perturbative series for quantum electrodynamics was argued on physical grounds by 



* sbravyi@us.ibm.com 

t divince@watson.ibm.com 

t daniel.loss@unibas.ch 

§ bterhal@gmail.com 



2 



Dyson [lOj]. There is a widespread belief that general perturbative series for quantum field theory and many-body 
physics do not converge but are to be viewed as asymptotic series, meaning that their lowest-order terms are a good 
approximation to the quantity of interest (while inclusion of higher-order terms may actually give a worse result, see 
e.g. i.) 

In this Letter we make these beliefs more rigorous by developing a formalism that justifies application of perturbation 
theory in the regime ||V|| ^> A. This formalism is applicable to many-body Hamiltonians that possess certain locality 
properties. A Hamiltonian H describing a system of n qubits will be called k-local if and only if it is represented as a 
sum of local interactions H — Hi such that each operator Hi acts on some subset of k or less qubits. In addition, 
the Hamiltonians in this paper will have the important property that each qubit occurs in at most a constant, say 
m = 0(1), number of terms Hi. The constraint says that for growing system size, the number of interactions in which 
any one qubit participates does not grow, but stays constant. Such a condition is generically fulfilled for Hamiltonians 
with shoft-range interactions studied in physics. An example is the standard (2-local) Heisenberg Hamiltonian on a 
lattice H — — £\ ■ J(rij)(XiXj + Y t Yj + ZiZj); the fact that the coupling J(r-y) is bounded-range, J(r^ > r) 0, 
ensures that each spin is acted upon by a constant number of terms. Note that our results are not restricted to 
Hamiltonians on lattices, but also hold for general networks of interactions such as expander graphs [ll[ . 

The interaction strength of a fc-local Hamiltonian H — Hi is defined as the largest norm of the local interactions, 
J = maxi || i?i || . Throughout this paper we assume that k is a constant independent of n. 

Our main result is a rigorous bound on the error that is made when a perturbative expansion for the effective 
Hamiltonian is truncated at a low order (the second or the third) in the regime where the interaction strength of V 
is small compared to A but ||V|| 3> A. This bound allows us to prove the following result. 

Theorem 1. Let H talgct be a k-local Hamiltonian acting on n qubits with interaction strength J. For any fixed 
precision e one can construct a 2-local simulator Hamiltonian H acting on 0(n) qubits with interaction strength O(J) 
such that the ground state energy of H approximates the ground state energy of -fftarget with an absolute error at most 
eJn. 

Let us comment on a significance of Theorem [TJ First of all, it eliminates the essential shortcoming of the earlier 
constructions P, 0, 0, S S 0| mentioned above, namely, the unphysical scaling of the interaction strength in the 
simulator Hamiltonian H. In our construction the interaction strength of H is bounded by a constant (depending 
only on k, e, and J) which makes it more physical and, in principle, implementable in a lab. The price we pay for 
this improvement is that we are able to reproduce the ground state energy only up to an extensive error eJn. For all 
realistic physical Hamiltonians the ground-state energy itself is generically proportional to n J, and thus the relative 
error can be made arbitrarily small. Secondly, estimating the ground-state energy of a fc-local Hamiltonian up to 
some sufficiently small constant relative error is known to be an NP-hard problem even for classical (diagonal in the 
|0), |1) basis) Hamiltonians. This hardness of approximation follows from the PCP theorem (PCP=probabilistically 
checkable proofs), see Chapter 29 of [ll[ for a review. Proving an analogous hardness of approximation result for 
quantum Hamiltonians is a widely anticipated development, a so called "quantum PCP theorem", see Section 5 in [I3 |. 
If true, it would imply that reproducing the ground-state energy up to some sufficiently small relative error (which 
is what our simulation achieves) is enough to capture the computational hardness of the target Hamiltonian. We 
also believe that the simulation results reported in this Letter will be relevant in actually proving such a quantum 
PCP theorem. Finally, the techniques we develop may be relevant in the context of adiabatic quantum computation 
where the simulator Hamiltonian must be capable of reproducing ground-state expectation values of local observables, 
see jjj. Il3j. 

The proof of Theorem 1 relics on a new way of deriving an effective Hamiltonian, which we call the Schrieffer- 
Wolff transformation after its use in many-body physics [l4(. Specifically, given an unperturbed Hamiltonian Ho, 
projectors P, Q onto the low-energy and the high-energy subspaces of Ho, and a perturbation V, the Schrieffer- Wolff 
transformation is a unitary operator e s with anti-hermitian S (S = —S^) such that (i) S is a block-off-diagonal 
operator, namely PSP = and QSQ — 0; (ii) the transformed Hamiltonian e s (Ho + V) e~ s is a block-diagonal 
operator, namely, Pe s (Ho + V)e~ s Q — 0. Given such an operator S, one defines the effective Hamiltonian on the 
low-energy subspace as H c s = P e s (Ho + V)e~ s P. A detailed exposition of the Schrieffer- Wolff formalism will appear 
in 0. 

Let us briefly sketch the proof of the theorem. The simulator Hamiltonian H is constructed using perturbation 
gadgets introduced in [JQ. "Gadget" is a technical term used broadly in theoretical computer science; in the present 
application, a gadget is simply an extra mediator qubit, and a Hamiltonian coupling the mediator qubit with some small 
subset of system qubits. For every mediator qubit u we define a projector onto its low-energy subspace P u = |0)(0| u 
and its high-energy subspace Q u = |l)(l| u . The purpose of a gadget is to simulate some particular fc-body interaction 
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-^target m the decomposition -fftarget = ^target + ^cisc- Here -ffeise are additional terms in the target Hamiltonian 
that we do not wish to treat using perturbation theory (since they are already 2-local, for example). 

Such simulation is achieved by applying perturbation theory to each gadget individually. A gadget's simulator 
Hamiltonian is 

H u = Hq+V u , where Hq = AQ U penalizes the mediator qubit for being in the state |1), and V u is a perturbation. 
With the proper choice of V u the effective Hamiltonian on the low-energy subspace, in which the mediator qubit u 
is in the state |0), approximates H^ aTget with an error e. Furthermore, this effective Hamiltonian can be obtained 
from H u via an approximate Schrieffer- Wolff transformation, that is, -ffj^ rgct — P u eS {Hq + V u ) e ~ S P u + 0(eJ) 
for some anti-hermitian S u satisfying p u S u P u = Q U S U Q U = 0. 

Combining the local gadgets together we get as a candidate for the simulator Hamiltonian H = Hq + V + -ffcisc with 
H n = J2 U Pq an d V — J2 U V u ■ Let A(-ff) and A(-fftarget) be the ground-state energy of H and -fftarget respectively. We 
prove that A(-ff) approximates A(iJ ta rgct) with a small extensive error by constructing a global unitary transformation 
e s mapping H to -fftarget (with a small extensive error), that is, -fftarget = P e s He~~ s P + O(enJ) where P = (g) u P u 
projects onto the subspace in which every mediator qubit is in the state |0). Given such a transformation one 
immediately gets an upper bound X(H) = X(e s He~ s ) < X(P e s He~ s P) — X(H talgc t) + 0(enJ). Here we have taken 
into account that restricting a Hamiltonian on a subspace can only increase its ground-state energy Making a 
natural choice S = 52 S u we prove that P e s He~ s P contains the desired term -fftarget and some cross-gadget terms 
where S u acts on Hq+V v with Using the 'independence' properties of the gadgets and the block-off-diagonality 

of S we are able to show that the contribution of these cross-gadget terms are small enough to be absorbed into the 
error term 0(ne), see Section II. 

In order to prove a matching lower bound, consider the transformed Hamiltonian H u = e Su H u e~ Su , where H u = 
Hq + V u . We prove an operator inequality H u > I u ® -ff t " rgct + O(eJ). Here I u is the identity operator acting on the 
mediator qubit u and O(eJ) stands for some operator with norm O(eJ). Intuitively one should expect this inequality 
to be true since the P-block of H u approximates H^ al . gct with an error O(eJ), the Q-block of H u contains a large 
energy penalty A, and the off-diagonal blocks PH U Q are small by the definition of the Schrieffer- Wolff transformation. 
Using the unitarity of e s and smallncss of S u we transform the above inequality into H u > 7" ® -fftarget + 0{eJ) 
which implies H > I <g> -fftarget + 0(neJ) and thus gives A(-ff) > A(-ff ta rget) + 0(neJ). These arguments are filled in 
at the beginning of Section II. 

Our results will be stated for two different mappings -fftarget — * H corresponding to two different gadgets, the one 
reducing the locality parameter k by a factor of 2, and the other reducing k = 3 to k — 2. By composing these 
mappings we arrive at Theorem [TJ 

Before proceeding with the details, let us make a few remarks about generalizations of this technique. We expect 
that similar results can be obtained for other perturbation gadgets in the literature (see e.g. [E Hi El ) j since these 
gadgets are all designed to work independently, i.e. a term in the target Hamiltonian is replaced by some local 
perturbed Hamiltonian H and cross-gadgets terms should have small contributions in the perturbation expansion. 
For Hamiltonians H of direct interest in many-body physics, we expect that it will also be possible to identify a gadget 
sub-structure H = J2 U H u sucn that H u gives rise to -ff t " rgot and cross-gadget contributions are small. For generic 
perturbed Hamiltonians, cross-gadget terms may not be small, which implies that the effect of the perturbations must 
be analyzed globally. We will consider such an analysis in a future paper 

Now we state two Lemmas used in the proof. The two Lemmas together can be regarded as an infinitesimal 
version of the Lieb-Robinson bound that governs time evolution of a local observable under a local Hamiltonian, see 



Lemma 1. Let S be an anti-hermitian operator. Define a superoperator L such that L(X) = [S,X] and L°[X] = X. 
For any operator H define r a (H) = \\e s He~ s \\ = \\H\\, ^i(-ff) = \\e s He~ s — H\\, and 



e.g. mm. 




p=0 



k > 2 



(1) 



where \\ ■ \\ is the operator norm. Then for all k > one has 



r k (H)<-\\L k (H)\\. 



(2) 



Using this Lemma one can show that 
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Lemma 2. Let S and H be any 0(l)-local operators acting on n qubits with an interaction strengths Js and Jh 
respectively. Let each qubit be acted upon non-trivially by 0(1) terms of H and S. Then for any k = 0(1) one has 

\\L k (H)\\=0(n-J k s J H ). (3) 

The proofs are rather elementary and can be found in the Appendix. 



I. THE GADGETS 

We shall use two types of gadgets proposed in H, namely, the subdivision gadget and the 3-to-2-local gadget. The 
former will be used to break fc-local interactions down to 3-local interactions while the latter breaks 3-local interactions 
up into 2-local interactions. A notational comment: in this section, all Hamiltonian operators, and Schricffer- Wolff 
operators S, refer to a single gadget u, and in the remainder of the paper these operators appear with the label u; in 
this section only, for economy, we omit this label. 

Subdivision Gadget. Let the target Hamiltonian be a single fc-body interaction H talgc t — JAB where A, B act 
on non-overlapping subsets of [fc/2] or less qubits and ||A||,||B|| < 1. Introduce one mediator qubit u, choose a 
parameter A> J, and define the simulator Hamiltonian H = Hq + V with 



ff = A|l)(l|„, V= ^AJ/2X U ® (-A + B) + V extIa . (4) 

Here T4 x t ra = {J/^){A 2 + B 2 ) acts trivially on the mediator qubit. Note that H contains only ([fc/2] + l)-body 
interactions. The purpose of the term V is to induce transitions \0} u — > \l) u — ► |0) M in the second order of perturbation 
theory, see Fig. 1, such that the corresponding effective Hamiltonian is proportional to (— ^4 + B) 2 containing the 
desired term AB and unwanted terms A 2 , B 2 which we cancel by V extr a- Next we define 

S = -i^ K Y u ®(-A + B). (5) 

One can check that S is the Schrieffer- Wolff transformation truncated at the first order, i.e. S — Si where Si is defined 
in Eq. (|23|) in the Appendix. We can calculate Pe He~ s P using the expansion in Lemma [T] A straightforward 
calculation shows that 



e s He~ s =(h+ [S, H] + ~{S, [S, H}]\ + 0( J 3 / 2 A' 1 / 2 ) = 



-^target 

AI+0(J) 



+ 0(j2A-5), ( 6 ) 



where we used that [S, T4 x t ra ] = 0. The upper and lower blocks correspond to the subspaces P and Q = I — P 
respectively. Thus, Pe s He~ s P is close to H tal gct — JAB, as desired; the error can be made 0{eJ) by choosing 
A = Je- 2 . 

3-to- 2-local Gadget. Let the target Hamiltonian be a single 3-body interaction, H targc t = J ABC, where A,B,C 
are one-qubit operators acting on different qubits and ||A||, ||£?||, ||0|| < 1. Introduce one mediator qubit u, choose 
A ^> J and define the simulator Hamiltonian H — Hq + V with Ho = A |l)(l| u , 

A - T- 

V = V d + V od + V cxtia , V d = -Aiji\l)(l\ u ®C, V od = -1^X U ®{-A + B), (7) 

and F ext r a = A 1 / 3 J 2 / 3 (-A + B) 2 /2 + J(A 2 + B 2 )C/2. Note that H contains only 2-body interactions. The purpose 
of the term V is to induce transitions |0)„ — > |1)„ — > |l) u — * |0)„ in the third order of perturbation theory, see Fig. 1, 
such that the corresponding contribution to the effective Hamiltonian is proportional to (— A + B) 2 C which coincides 
with ABC up to some unwanted terms which are canceled by V ex t Ta . We define 



iX 

S=—Y u ®{-A + B) 



I + xC + x 2 C 2 -^(-A + B) 2 



(8) 



One can check that S is the Schrieffer- Wolff transformation truncated at the third order, i.e., S — Si + S2 + S3, see 
Eq. (|23p in the Appendix. We calculate the effective Hamiltonian Pe s He~ s P. Let us first estimate an error resulting 
from cutting off the expansion, see Lemma [TJ Recalling that L — [S, ■] one gets 

\\Pe s He- s P -P(H+ [S, H] + \[S, [S, H}]\ P\\ < \\\P L%H) P\\ + ^\\L\H)\\. (9) 
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l7> 



lo> lo> 

FIG. 1: Transitions induced by the perturbation V for the subdivision gadget (left) and the 3-to-2-local gadget (right). 



We note that P L 3 (H) P = P L 3 (V d) P since each application of S flips the mediator qubit and a non-zero contribution 
comes only from the terms with an even number of flips. Using a bound \\S\\ = O(x), see Eq. ([8]), one can upper-bound 
the r.h.s. in Eq. (|9|) as 0( J 4 / 3 A -1 / 3 ). A direct but lengthy calculation shows that 



e b He 



H + [S, H] + i[S, [S, H}] + [S, [S, H]}]\ + O^' 3 ^ 1 ' 3 ) 



^target 

A/ + 0(A 2 / 3 J 1 /3) 



where H*. 



. (j4/3 A -l/3 )) (10) 

J ABC, as desired. The error 0(J 4 / 3 A -1 / 3 ) in H targot can be made 0(eJ) by choosing A = J(T 3 . 



II. COMBINING THE GADGETS TOGETHER 



Let i?tar g ct = J2u ^target +-^elsej where each term ii t " rgct can be dealt with using one of the gadgets described above. 
The simulator Hamiltonian is H = H u +H e \ se where H u is the simulator constructed for a gadget u as above. Using 
Eqs. (|5irU|) one gets the inequality e s " if e" s " > I u ®H^ Igct +0(eJ) which yields H u > e~ s " (I u <g>H% aT&et )e su +0(e J). 
Applying Lemmafflone gets ||e- s "(f g) if" arget )e s " - f ® H^ Igct \\ < \\ [S«,f <8 ff« rget ] || = O(eJ), see Eqs. ((511). 
It follows that if > f (gi flta, r g e t + 0(eJ). Summing up these inequalities over all gadgets one arrives at i? > 
7 (g) -fftarget + O(enJ), where I acts on the mediator qubits. Thus 

\{H) > X(H target ) + O(enJ). (II) 

Bounding Cross-Gadget Contributions. Let P = §§ u P u be the projector on the subspace in which every 
mediator qubit is in the state |0). Define S = J2 U ^ u where S u is constructed using Eqs. (|5l8p . Below we shall prove 
that 



\\Pe s He~ s P-H targct \\=0{enJ). (12) 

Then one can get an upper bound on X(H) by observing that 

X(H) = X{e s H e' s ) < X(P e s H e~ s P) = X(H target ) + 0{enJ). (13) 

Combining Eqs. (|11I13|) one gets \X(H) — A(i?target)| = 0{enJ) which is the desired result. It remains to prove 
Eq. fT2")) . The idea is to bound the cross-gadget terms using Lemma [1] and Lemma [5] Two important properties (valid 
for both gadgets) will be used repeatedly in this argument: (1) S u always flips a mediator qubit it; (2) For u ^ v one 
has [P u , if] = 0, [P u , S v ] — 0, [Hq, S v ] — 0. This latter property essentially captures the independent action of the 
local gadgets. 

Let us start with the subdivision gadget. Recall that we have chosen A = Je~ 2 . From Eq. we can see that S is 
a 0(l)-local operator with an interaction strength 0(e). Applying Lemma [T] we get 

\\Pe s He- s P- i7 taxget || < \\ [S,H else ] || +£||P(if + [S,if ] + ±[S, [S, if]]) P - H^. gct \\ + ^||L 3 (fl)||, (14) 

u 

Lemma [5] gives 1 1 [S, H e i se ] | 1 = O(enJ) and ||L 3 (i7)|| = 0(e 3 An) = O(enJ). Properties (1),(2) above imply 
P[S V ,H U ]P = for u ± v and that P [S Ul , [S" 2 , if]] P = unless Ui = u 2 . Property (2) implies that 
P [S v , [S V ,H%]] P = for u ^ v and thus the only non-zero cross-talk term could be P [S v , [S v , V u }} P where u ^ v. Us- 
ing the explicit form of V u , see Eq. @, and property (1) one concludes that P [S v , [S V ,V U ]] P = P [S v , [S v , F e ^ tr J] P. 
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By construction of the individual gadgets, the contribution of the "diagonal" terms (those in which H u and S u belong 
to the same gadget) is O(enJ), see Eq. ©. Combining these observations together we arrive at 

\\Pe s He- s P - H taTset \\ < 0(enJ) + £ || P [S\ [S\ V^]] P || (15) 

Now, note that for any u there exist only 0(1) mediator qubits v such that [S v , Kxtra] 7^ 0. S ince V^tra has norm 
O(J) the contribution of the cross-talk terms is 0(ne 2 J) and we arrive at Eq. (|12p . 

Let us do a similar analysis for the 3-to-2-local gadget. Recall that we have chosen A = Je~ 3 . From Eq. (|HJ) we see 
that S is a 0(l)-local operator with interaction strength 0(e). Applying Lemma[T]we get 

|| Pe S He- S P - fftargetH < || [S, H clsc \\+J2\\P + [S, H U ] + \[S, [S, H U ]}) P - target II 

U 

+ L\\PL 3 (H)P\\ + ±\\L\H)\\. (16) 

Lemma |U gives || [S, H e i se \\ = O(enJ), \\L 4 (H)\\ = 0{e 4 An) = O(enJ). Property (1) implies PL 3 {H )P = 0. 
Taking into account that V has interaction strength A 2 / 3 J 1 / 3 = Ae, see Eq. ([7]), and applying Lemma H we get 
\\PL 3 (H)P\\ = \\PL 3 (V)P\\ = 0(ne 3 Ae) = O(enJ). Repeating the same arguments as for the subdivision gadget we 
conclude that the only cross-talk terms contributing to Eq. (fT6|) are P [S v , [S v , V u ]] P where «/». Using the explicit 
form of V u , see Eq. ©, one concludes that P [S v , [S v , V u ]] P = P [S v , [S v , V^ tra ]j P so we again arrive at the bound 
Eq. (p~5|) . By definition, T4" tra has norm 0(A 1 / 3 J 2 / 3 ) = (9(Je _1 ) and S v has norm 0(e). Making use of the properties 
of i/target we conclude that the contribution of the cross-talk terms is O(neJ) and we arrive at Eq. (fT2"]) . 

III. SIMULATION OVERHEAD 

Let us estimate the overall increase of the interaction strength associated with m levels of simulation using the 
subdivision gadget. Let Hi be a Hamiltonian at a level i of the simulation such that Hq = H ta , Tge t and H m is 3-local. 
Accordingly, Htarget can be A:-local where k ~ 2 m . The overall number of mediator qubits one needs for the simulation 
is 0(n2 O( - fe - ) ). For simplicity one can add all these qubits to the system from the beginning by letting H ta , rgc t to act 
trivially on them. It allows us to assume that at every level of simulation we have n' — 0(n2 ^) qubits. Let A.; 
be the ground state energy of H^. Let Ji be the interaction strength of Hi. In order to make |Aj+i — A^| ~ n'S it 
suffices to choose J l+ ± = 8~ 2 J 3 . Accordingly, S^ 1 J rn = Jo) 3 ™ and thus J m =o(io g k) ~ Jq° 1v ^ S~ poly ^ . Choosing 
S = 2~ u[ - >e one can make the overall error |Ao — A m | ~ en. Assuming that J = 0(1) we get J &nal ~ (2e- 1 Y° ly{ ^ . 

It may be possible to improve this exponential scaling by using a direct fc-to-2-local gadget. 
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APPENDIX 

Proof of Lemma\^ For any real t > 0, let H(t) = e St H e~ st . Define an auxiliary quantity 

r k (H,t) = \\H(t)-J2^LP(H)\\. (17) 

Let us get an upper bound on the increment r% (H, t + St) — r k (H, t) . First of all notice that 

H(t + St) - H(t) = St e st L(H) er st + 0((St) 2 ). (18) 
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Secondly, for any k > 1 one has 

g (H-fl)* ^ _ g ^ LP(iJ) = 5t g ^ (i(fl)) + 0((ft)2) (lg) 

p=0 p=0 " p=0 

(if fe = 1 then there are no 0(St) terms on the right hand side). Applying the triangle inequality one gets 

r k (H,t + St)-r k {H,t) < Str k ^(L{H),t) + 0((St) 2 ). (20) 

Since r k (H, i) is a continuous function of t, it is legitimate to add up the inequalities Eq. (|20p for t = 0, St, 2St, . . . , s, 
take a limit St — > 0, and replace the resulting sum by an integral, which yields 

r k (H,s)< f dtr k -i(L(H),t) if k > 1, and r {H, s) = \\H ||. (21) 
Jo 

Applying this upper bound recursively and evaluating the integrals one arrives at 

t k 

r k (H,t) < \\L k {H)\\ ■ — forallfc>0. (22) 
(k)\ 

Since r k (H) = r k (H, 1) Eq. ([22]) proves the lemma. □ 

Proof of Lemma\^ Represent S and H as a sum of local operators S — X)j=i &i an d H = ^2f=i such that any Si, 
Hj act on 0(1) qubits and any qubit is acted on by 0(1) operators Si, Hj. Let us use the term elementary commutator 
of order k for a multiple commutator that involves some Hj and k operators Si. For example, [Si, [S^ifi]] is an 
elementary commutator of order 2. Note that a commutator of any 0(l)-local operators is again a 0(l)-local operator. 
Therefore for any constant k the number of non-zero elementary commutators of order k contributing to L k (H) can 
be bounded as 0(n). Each elementary commutator of order k has a norm at most 2 fc j| Jh- Applying the triangle 
inequality to \\L k (H)\\ we get Eq. ©. □ 



SYSTEMATIC SOLUTION FOR THE SCHRIEFFER- WOLFF TRANSFORMATION 

One can note that the equation P e s (Ho + (3 V)e~ s Q — has a unique solution in terms of a formal Taylor series 
S = YlpLi Sp P P 'i tne series coefficients S p can be straightforwardly computed. In the gadgets we have used an 
approximate version of S obtained by truncating the Taylor series at the first (subdivision gadget) or the third order 
(3-to-2-local gadget). Introducing operators Vd = PVP + QVQ, V d = PVQ + QVP and linear maps L = [Ho,-], 
L p = [S p , ■] one can derive that 

Si = L ~ 1 (V o d), S 2 =L- 1 L 1 (V d ), S^-^L^LUH^+L^^Vd). (23) 



[1] J. Kempe, A. Kitaev, and O. Regev. The Complexity of the Local Hamiltonian Problem. SIAM Journal of Computing, 

35(5):1070-1097, 2006. Earlier version in Proc. of 24th FSTTCS. 
[2] D. Aharonov, W. van Dam, J. Kempe, Z. Landau, S. Lloyd, and Oded Regev, Adiabatic Quantum Computation is 

Equivalent to Standard Quantum Computation SIAM Journal of Computing, 37(1):166-194, 2007. 
[3] R.I. Oliveira and B.M. Terhal. The complexity of quantum spin systems on a two-dimensional square lattice. 2005, 

http : // arxiv . org/abs/quant-ph/0504050 
[4] S. Bravyi, D.P. DiVincenzo, R. Oliveira, and B.M. Terhal. The complexity of stoquastic local Hamiltonian problems. To 

appear in Quant. Inf. and Comp., Vol. 8, No. 5, p. 361, 2008. 
[5] S. Jordan and E. Farhi. Perturbative gadgets at arbitrary orders, http://arxiv.org/abs/0802.1874 
[6] Norbert Schuch and Frank Verstraete. Interacting electrons, Density Functional Theory, and Quantum Merlin Arthur. 

|http : //arxiv . org/abs/0712 . 0483 | 
[7] Jacob D. Biamonte and Peter J. Love. Realizable Hamiltonians for Universal Adiabatic Quantum Computers 

|http : //arxiv. org/ abs/0704. 1287 | 
[8] M. Reed and B. Simon. Methods of Modern Mathematical Physics IV: Analysis of Operators. Academic Press, San Diego, 

CA,USA, 1978. 



8 



[9] S. Bravyi, D. DiVincenzo, and D. Loss, Polynomial-time algorithm for simulation of weakly interacting quantum spin 
systems, http : // arxiv . org/abs/0707 . 1894 
[10] F. J. Dyson. Divergence of perturbation theory in quantum electrodynamics. Phys. Rev., 85(4):631-632, Feb 1952. 
[11] V. Vazirani Approximation Algorithms Springer 2001. 

[12] D. Aharonov, D. Gottesman, S. Irani, and J. Kempe The power of quantum systems on a line. Proc. 48th FOCS, p. 

373-383 (2007). |http : //arxiv . org/abs/0705 ■ 4077] 
[13] Stephen P. Jordan, Edward Farhi, and Peter W. Shor. Error-correcting codes for adiabatic quantum computation. Phys. 

Rev. A, 74:052322, 2006. 

[14] J. R. Schrieffer and P. A. Wolff Relation between the Anderson and Kondo Hamiltonians, Phys. Rev. 149 p. 491, 1966. 
[15] A. Kitaev. Anyons in an exactly solved model and beyond. Annals of Physics, 321:2, 2006. 

[16] S. Bravyi, D. P. DiVincenzo, D. Loss, and B. M. Terhal. Perturbation theory of many-body systems and the Schrieffer-Wolff 

transformation. In preparation. 
[17] M. B. Hastings and T. Koma. Spectral gap and exponential decay of correlations. Comm. Math. Phys., 265:781, 2006. 
[18] S. Bravyi, M. B. Hastings, and F. Verstraete. Lieb-Robinson bounds and the generation of correlations and topological 

quantum order. Phys. Rev. Lett., 97:050401, 2006. 
[19] By abuse of notation, P e s He~ s P denotes the restriction of the operator e s He~ s onto the subspace P. 



